I 



X 



I. INTRODUCTION 



Antiferromagnetism and Superconductivity in UPt 3 

Anupam Garg 

Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208 

(today) 

OO 

On 

' The short ranged antiferomagnetism recently seen in UPt3 is proved incompatible with two di- 

mensional (2D) order parameter models that take the antiferromagnetism as a symmetry breaking 
» I . field. To adjust to the local moment direction, the order parameter twists over very long length 

scales as per the Imry-Ma argument. A variational solution to the Ginzburg-Landau equations is 
used to study the nature of the short range order. Although there are still two transitions, the lower 
one is of first order — in contradiction to experiments. It is shown that the latent heat predicted by 
l/") ' the 2D models at the lower transition is too large not to have been seen. A simple periodic model 

is numerically studied to show that the second transition can not be a crossover either. 
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The purpose of this paper is to discuss the implications of recent experiments (lj on antiferromagnetic (AFM) order 
in UPt3 for the superconducting state. It is found that (i) a very weak AFM that sets in at 5 K is multi-domain in 
structure with a domain size of about 300A, and (ii) despite a magnetic moment per U atom of only 0.02 /xg, the 
domains can neither be moved nor oriented by magnetic fields upto 3.2 T. We shall show that these findings pose 
'-^J \ great difficulty for certain widely studied models in which the superconducting order parameter transforms as 
£5 ■ one of the two dimensional (2D) represenations of the relevant point group, D' 6h , and in which the AFM (or even 
less plausibly, a structural modulation) provides a symmetry breaking held that splits the superconducting transition 
, ^ , , into two just as an applied magnetic field splits the A transition in superfluid 3 He. We believe that these data make 
it overwhelmingly likely that superconducting UPt3 is described by two order parameters belonging to independent 
representations HH, although the precise representations involved are still debatable. 

Two main experimental facts with which any theory of UPt3 must agree are: (A) There are two superconducting 
transitions at T c+ and T c _ in zero field T c+ ~ 550 mK, and AT C = T c+ — T c _ ~ 50 mK. Both transitions appear 
to be of second order, with no observable latent heat. (B) For basal-plane magnetic fields, H 1 c, the phase diagram 
£^ \ is isotropic 0; i.e., it is the same for all directions of field.Q We deliberately do not exploit other facts that also bear 
. on the order parameter symmetry in this paper. Chief among these is that for H 1 c, there is a tetracritical point in 
the magnetic field-temperature (H-T) plane where four second order lines meet ||]. The original 2D models g| are 
0^ i incompatible with a tetracritical point unless H 1 c |,|l^|, which also argues against them, as the observed phase 
diagrams for H||c and H at 45° to c are similar to those for HI c P,|lT|, although the tetracritical points in these 
cases may only be apparent ones due to limited resolution of the phase boundaries ||[l2| . Secondly, a specific E2 U 
model is compatible with a tetracritical point for all field orientations |3| . It is therefore useful to have independent 
arguments that constrain the order parameter without relying on the presence or absence of tetracritical points. 

With this motivation, we analyze in this paper the Ginzburg-Landau (GL) theory for 2D models jHH] when the AFM 
symmetry breaking field is allowed to vary randomly in direction over a length scale a = 300A.gSince the problem 
involves an order parameter with a continuous symmetry in the presence of a random field, the Imry-Ma argument G3] 
implies that the long-ranged average of the order parameter itself must vanish. Other forms of order are not ruled 



1 We shall ignore an anisotropy of 4-5% H as it is too small to affect our arguments. 

2 To avoid misunderstanding, we stress that we are not studying the AFM transition itself, or the reason for AFM disorder 
Jl3[ . Since the AFM order sets in well above the superconducting transition, it is justified to treat it as an external field. 
The spatial variation of this field in no way invalidates the usual argument for using GL theory to interpret superconducting 
transitions. Briefly, the argument (also known as the Ginzburg criterion) is that the critical region is of width (d/£o) relative 
to To, where d is an inter-electron distance, and £o is the coherence length. For UPt3, d ~ 5A, and £o ~ 200A, so this region 
is unobservably small. In particular, the question whether the transition is sharp or smeared out by the random AFM is 
experimentally irrelevant. 
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out, however, and a balance of gradient and field orientation energies shows jl5| that there are two transitions, the 
upper one of second order, and the lower one of first order. If one is only interested in global properties of the phase 
diagram and the universality class of the transitions, there is nothing more to be said. We are more interested in 
whether or not this scenario could describe the experimental data on UPt3. To this end we calculate the latent heat 
at the lower transition assuming that the coupling to antiferromagnetism is weak. We find that this is too large to 
have been missed by the experiments |(| if the splitting is taken to equal the observed value. In other words, we can 
dismiss the possibility that the lower transition seen experimentally is of first order with an unobservably small latent 
heat. A second possibile explanation for the data (which arises when the AFM coupling is not weak) is that the two 
transitions predicted by the theory are too close to be resolved, and that the lower specific heat jump is actually due 
to a cross-over and not a true transition. We exclude this possibility too. Using a simple periodic model for the AFM 
domain structure, we find that the observed H C 2 slopes ||[llj and the 300 A domain size imply that a sufficiently 
sharp cross-over must be accompanied with a AT C of order 80 K, which is clearly absurd. 

To orient the discussion further, let us review how the 2D models account for facts A and B when the AFM is 
uniform. The GL free energy /gl for all these models can be written in terms of a vector if — (rj Xl rj y ) with complex 
components. In zero magnetic field /gl is the sum of bulk, gradient,^ and symmetry breaking field terms: 

/bulk = —(T - T ) v * Vl + + fclmrhl 2 , (1-1) 

/grad = Kxditfdirij + K 2 d t r]*djr]j + Kzditfdj^ 

+K 4 d z T)*d z r)i, (1.2) 

/ S BF-2e(|n^| 2 -, / *^/2). (1.3) 

Here, i,j € {x,y}, and the summation convention is used, n is the basal-plane unit vector along the magnetic 
moments, and e, (3i, and /?2 are all positive. (If 02 < 0, there is only one transition, and the model can not describe 
UPt3.) Suppose n = y everywhere. Then the first transition occurs at T c+ = Tq + e to a phase with fj = e^(l,0). 
Ignoring the global phase cf>, we refer to this as the real or (1,0) phase. It breaks the rotational symmetry of the 
normal state. The second transition occurs at T c _ = To — /3ie//?2 to a phase with fj — e l( ^(l,iu) where u is real and 
grows smoothly from to ±1 as T is lowered further. This phase breaks rotational and time reversal symmetry, 
and we shall refer to it as the axial or the (1, i) phase. This explains fact A. To explain fact B, it is necessary in 
the original 2D models Q to assume that the symmetry breaking field rotates in response to the applied magnetic 
field. A symmetry breaking perturbation that was locked to the lattice would split the T c 's, but would not yield an 
isotropic basal-plane phase diagram. In particular, the tetracritical point would not exist for some H orientations. 
This assumption was justified by hypothesizing [Jl5| that the AFM anisotropy energy was so weak that the Neel vector 
reoriented itself to stay perpendicular to the magnetic field even for relatively small H, and this idea has even been 
refined (li| to explain the minute anisotropy of H C 2 found by Keller et al. ||. 

The neutron scattering data of Lussier et al. jl| directly refute the above hypothesis. They find that an AFM 
domain is of average size a ~ 300 A ^ and is associated with one of three choices for the spin density wavevector q, 
with n||q. The scattering intensity under the peaks with different q's is essentially unchanged by fields of up to 3.2 
T, for either H||q or H _L q, and in both field-cooled and zero-field-cooled experiments. In addition, application of 
a field causes no transfer of intensity to a scattering wave vector Q which is forbidden (because Q x n = 0) when 
H = 0. This shows that the moments do not cant away from the q direction while preserving the wavevectors. 

Why, given these data, should one further investigate the 2D models? One reason is the modified 2D model H . In 
this model, k 2 = «3 = 0, and the AFM is assumed to split the K\ term into n^\dir^ | 2 + Ki\diT]±\ 2 , where ryy and r]± 
are components of fj parallel and perpendicular to n, and — n7 = 0(e). (A similar split ting is assumed for K4.) 
Note that the AFM continues to couple directly to the order parameter as described by Eq. ( |l.3| ). The linearized GL 
equations for rj\\ and rj± then decouple even when H ^ 0, and an isotropic basal-plane phase diagram is recovered 
even if the AFM is rigidly locked to the lattice. However, the random AFM excludes this model too, as the scaling of 
gradient and field orientation energies with e, and the overall balance between them, is unchanged. A second reason 
is that although the general thrust of the Imry-Ma argument [jl4| must continue to hold, the problem for UPt3 is 
not directly equivalent to an rt-component vector spin model, and we are unaware of any specific analysis for it. In 



3 The gradient energy in the modified 2D model |^] is slightly different, as discussed below. 

4 It follows that a structural modulation which should, a priori, be insensitive to the orientation of H, also can not be a viable 
symmetry breaking field. 

5 The absence of anomalies in heat capacity Jl^ or NMR linewidth, Knight shift, or 1/Ti Jfij], raises doubts whether the 
AFM is static long range order; we set these doubts aside here as they make the 2D models even more problematic. 
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particular we know of no careful comparison of the energetics of the axial phase versus a real state with a wandering 
axis of orientation. 

The rest of the paper is organized as follows. The GL free energy is analyzed in Sec. II. We perform a variational 
analysis to understand the nature of the short-range order at the upper transition, and also calculate the latent heat 
at the lower transition. The simple periodic model is analysed in Sec. III. A brief discussion in Sec. IV concludes the 
paper. 



II. ANALYSIS OF THE GL FREE ENERGY 



We wish to minimize the GL free energy, which is still given by Eqs. (1.1-1.3), but n(r) is now r dependent with a 
correlation length a — 300 A. We first simplify the gradient energy to 



/grad = K\Vr]i 



(2.1) 



as we are interested only in the competition between the gradient and the symmetry breaking energies, and not in 
the precise slopes of the H C 2 lines and whether they cross or not.^ It is convenient to write n = (cos 0, sin ff), and 
introduce the gradient energy scale Vd = n/a 2 , and circular polarization components, T]l,r — 2~ 1 ^ 2 (ilx ^ 

The first step is to find the critical temperature T c , and the nature of the critical order parameter. To do this we 
drop the quartic terms in /gl, and minimize the rest. This yields an eigenvalue problem given by the linearized GL 
equations: 



ee™ m 



-tVR, 



(2.2) 
(2.3) 



where t = T c — Tq is the T c enhancement. The critical so luti on is t hat with the largest value of t. 

Two qualitatively different types of solutions to Eqs. (2.2) and (2.3) are as follows. The first, or local solution, is 



motivated by the answer for uniform 6 . Then ff is also uniform, the gradient terms vanish, and we get t]l — —e 2l0 r]ii, 
t = e. If we now require ff _L n, i.e., r/^ = —e 2lS r]ii, even when 9 is varying, then we pick up a gradient energy of the 
order of Vd, and ii OC ai = e — ctVd, where a = 2a 2 ((X78) 2 ) is a number of order one, and (•) denotes a spatial average. 
This solution does well only when e ^> Vd, i.e., when the symmetry breaking energy is dominant. The second, or 
real solution, has the form ff = e^m, m being a real unit vector varying on a length L»n. The gradient energy is 
clearly of order Vd(a/L) 2 . A random walk argument for symmetry breaking energy leads to the estimate — e(a/L) 3//2 . 
Minimizing the sum of these two energies, we obtain L ~ a(Vd/e) 2 , and t rca i ~ e(e/Vd) 3 . This solution requires Vd 3> e 
for consistency, and the T c enhancement is strongly reduced compared to the uniform case. 

A correct estimate of the symmetry breaking energy for the real solution is more subtle, however. Although long 
range order in ff is forbidden, it is not obvious that the estimate of the energy itself is correct, as one could try and 
lower it by short range adjustment to the local value of n. Indeed, we will find that though the T c splitting is governed 
by the e 4 /V d 3 energy gain, this gain is subdominant to that due to the short-ranged order, and the enhancement in 
T c itself is larger, of order e 2 /Vd- 

The above results follow from a variational argument, which we present in three stages. In the first stage, let us 
consider the trial wave function -q^ = 1, r]R = (e~ 2l6 ( r \ with £ as a variational parameter. We call this the axial 
solution. (See reasons below.) The optimal value of £ is to be found by minimizing 



W 



(v\Hgl\v) 

(v\v) 



(2.4) 



where Hql is the linear operator, or GL "Hamiltonian" in Eq 
obvious quantum mechanical notation for the averages. A simp 



|) and (2.3), is the trial solution, and we use an 

e calculation gives 



W 



: Re(Q + a^|C| s 
1 + ICI 2 



(2.5) 



6 This simplification is equally justified for the modified 2D model [jj. The gradient energy cost for variations of ff along and 
perpendicular to n is now slightly different, but it is, nevertheless, always positive. Any suitable average for the gradient energy 
scale Vd therefore suffices for the scaling arguments that follow. 
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[The quantity a = 2a 2 ((V 0) 2 ) as before.] Minimization gives a T c enhancement i a xiai = — ocVd + (e 2 +a 2 V d 2 ) 1 / 2 , which 
is w e 2 /2aVd for e <C Vd, and ~ e — aVd when e 3> Vd- We have thus done better than both the local and the real 
solutions. 

In the second stage, we generalize the above solution by taking 

VL = 1, m (r) = J d 3 v'Q{v - v')e- 2 ^"> = h{v), (2.6) 

where Q(r) is a variational kernel. We take e 2 " 10 ^ to be Gaussian-distributed, with [e 2%e ^) = 0, and a correlation 
function g(x) = ( e 2ie(r) e -2i0(r+x)^ decaying on a length scale a. Denoting the volume of the system by fi, the average 
of the symmetry breaking energy is found to be 

2 Re^ J d 3 r vlVRe 2ie(r) = 2eF 1 [Q], (2.7) 

where Pi[Q] is the functional 



d 3 k 
(2nj~ 3 



F ilQ}= I 7^3<3( k Mk), (2.8) 



where Q(k) and <?(k) are Fourier transforms of Q(r) and g(r). Similarly, the average of the gradient energy equals 

^ J d 3 r\V VR \ 2 = kF 2 [Q], (2.9) 



where 

d 3 k 



F 2 [Q] = I j^k 2 Q 2 (k)g(k). (2.10) 



Finally, the spatial average of |?7i?| 2 , needed for normalization, is given by 

d 3 k 



1 J d 3 r\ m \ 2 = F 3 [Q] = J pQ 2 (k) S (k). (2.11) 



nj 1 JL " J ~ y (2tt) 3 

The expectation of i?GL is therefore given by 

W=(2eF 1 [Q] + K F 2 [Q])/(l + F 3 [Q]), (2.12) 

which must be minimized with respect to Q(k). Setting the variation of W to zero, we obtain 

(1 + F 3 )(2eSF 1 + k5F 2 ) - (2eF 1 + kF 2 )SF 3 = 0. (2.13) 

Using 5Fi/5Q(k) = <?(k), 2k 2 Q{k)g{k) : and 2Q(k)g(k) for i = 1, 2, and 3, respectively, and cancelling g(k) in this 
equation, we get 

(1 + F 3 )[e + Kk 2 Q(k)} - (2eF 1 + K F 2 )Q(k) = 0. (2.14) 
This yields the Yukawa-like form 

where k 2 , is self-consistently given by 

In the last form, which follows by comparison with Eq. ( 2.12| ), W m i n denot es the minimum value of W. 

It remains to find the Fi[Q] with the form ( [2.15 ) and substitute in Eq. (2.16) to obtain fc and Wmin. Since g(x) 
decays on a length scale a, we expect g(k) to decay on a scale a -1 . Since g(x = 0) = 1, this decay must be faster 
than 1/fc 3 as k — ► oo; this fact will guarantee convergence of the integrals below. For F\ we obtain 
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e 1 



Fi = — - 



k2tt 2 



fc 2 + fc 2 



g(k)k 2 dk. 



(2.17) 



Now, by the second equality in Eq. (2.16), we expect ko to vanish as e — ► 0. It is therefore valid to put ko to zero in 
the integrand above. This yields 



Fx 



where 



1 



Cl 



g(k)dk 



2n 2 a 2 

is a dimensionless number of order unity. In the same way, we get 

k 2 



Fn 



K 2 2tt 2 J q (k 2 + k 2 ) 2 



g{k)k dk « Ci(ea/«) , 



where the last result follows from again setting ko = in the integrand. Lastly, F$ is given by 

1 



F* = — - 



k 2 2it 2 J Q (k 2 + k 2 ) 



-g{k)k z dk. 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



It is impermissible to set ko = now as the resulting integral is divergent. Instead, we note that the integrand is 
sharply peaked near k = 0, and use this fact to replace g(k) by go = g(k — 0). This yields 



F 3 « c 2 



1 

k a \Vd 



where c 2 — go/ 87ra 3 is another dimensionless number of order unity. 

We will show shortly that F 3 -C 1. (T his is also to be expected from the fact that F3 = (\t]r, 
be the case, using Eqs. ( [2.18| ) and ( 2.20 ), and neglecting F 3 in comparison with unity in Eq. ( 2.16 ), we obtain 



(2.22) 

2 ^ Assuming this to 



k Q a = v / cT(e/V r d ), 



(2.23) 



iaxiai = -W min = Cl e 2 /V d . (2.24) 

It then follows that F 3 ~ e/V d <C 1 self-consistently. The assumption fc a <§; 1 is also seen to be true. Although the 
T c enhancement is still of order e 2 /Vd, we can see that it exceeds that obtained with the earlier simple trial solution 
by noting that the quantity 2ac\ can be written as the product (k 2 )(k~ 2 ), where (•) is a mean defined with respect 
to the distribution g(k)d 3 fc/(27r) 3 . The desired result then follows from noting that 2ac\ > 1, i.e., c\ > l/2a, either 
as a consequence of Holder's inequality, or by recalling that the arithmetic mean of positive numbers with a non-zero 
spread always exceeds the harmonic mean. Further, the solution is non-perturbative since \rjn\ ~ (e/Vd) 1//2 , which is 
not an integer power of e. 

The improvement over the real solutio n is illusor y, ho wever. First, note that the axial solution is degenerate, since 



if (t]l,t1r) = i u i v ) is a solution to Eqs. (2.2) and (2.3), so is (v*,u*). (This is merely complex conjugation.) Since 
T]r 7^ rf L for Eq. ( |2.6| ), these are distinct solutions. On the other hand, the GL equations are seen to be completely 
real when written in terms of r\ x and rj y> suggesting that the ground state is also real and unique.^ Thus a linear 
combination of the solution just found with its complex conjugate might be even lower in energy. The second point 
is that the length scale fc^" 1 ~ e _1 is much less than the Imry-Ma length scale ~ e~ 2 . Thus the axial order found 
above is only short-ranged, and we can multiply the variational answer by a phase factor 7(1") varying on a length 
scale greater than k^ 1 and try to further reduce the energy in this way. This forms the third stage of our argument. 



We denote the solution (2.6) by fj^\ its complex conjugate by fj^ 2 \ and examine the combination 



7 Another way of saying this is that in a large but finite system, Hgl has no symmetries (including time-reversal) with a 
random AFM. Thus this argument does not apply to the periodic model studied in in Sec. III. 
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u(r) 



e -iT(r)^(2) 



(2.25) 



This has cartesian components 



— sin 7(1 + Re ft-) — cos 7 Im /i, 
cos 7 (1 — Re h) + sin 7 Im h. 



(2.26) 
(2.27) 



Since |/i(r)| <§C 1, this solution is close to (— sin7, cos 7), i.e., to a rotation of the ( 1, 0) solution. If 7 were a constant, 
the energy would be identical to that of the solution (2.6), i.e., given by Eq. ( 2.24 ). To calculate the difference, let us 
define 



Tl — Tq + taxial , 



(2.28) 



and A 2 = 1 + (\h(r)\ 2 ) = 1 + Fa. If we average the quadratic part of Jgl over a length scale of order k , we obtain 
a coarse grained free energy 



JC' 



A 



((T —T) + K (V 7 ) 2 - e (cos2[0(r) - 7 (r)]> Cfl ) 



(2.29) 



where 



denotes a coarse grained average. Note that the fi rst te rm is proportional to 7\ — T and not T — T as in 



Eq. (1.1). We can now apply the Imry-Ma argument to Eq. ( 2.2S ). Suppose the energy is minimized when 7 varies 
on a length scale L. In a supercell of linear size L, there are approximately (L/a) 3 cells in each of which 9(r) can be 
regarded as constant and independent of the other cells. Thus the coarse grained average of either cos 28 or sin26 l in 
a supercell is a number of order (L/a) 3 ^ 2 of either sign. By choosing cos 27 ex +(cos2#) cff , and sin 27 cx +{sm29) cg , 
we obtain a symmetry breaking energy which is systematically negative in every supercell. Since the gradient energy 
is of order L~ 2 , we reco ver th e Imry-Ma length scale L ~ e~ 2 , an d a n energy lowering of 0(e 4 ). This energy is 
subdominant to the gain ( 2.24 ) from the short ranged order of Eq. (2.6), but since it represents an additional gain, 
the true solution is indeed disordered and real. It is in this sense that the real solution must be understood. 
Therefore, the transition temperature for the real solution is given by 



T c+ = T 1+ 0(e 4 /V 3 ). 



(2.30) 



(Actually, T\ — Tq could also contain a term of order e 4 , but that is immaterial.) Let us now consider what happens 
at lower temperatures. Now the quartic terms come into play, and the positive /3% term in /gl favors an axial 
solution. For e < Vj, we can compare the free energies of the real and axial solutions directly. Since (|r?i?7i| 2 ) vanishes 
for a purely axial solution and equals ((rjiVi) 2 ) f° r a purely real one, these free energies will differ strongly in the 
contribution of the term. Approximating multiplicative factors such as (M 4 )/(M 2 ) 2 by unity, and denoting real 
and axial solution quantitites by subscripts r and a, the standard analysis of a Landau free energy with a second 
order transition yields 



f r = -C r (T- 

fa = -C a (T- 



- T c+ ) 2 /2T c+1 
-Ti) 2 /2Ti, 



where 



C r /T c+ = l/2(/3i +/3 2 ), 
C a /T x = l/2/3i . 



(2.31) 
(2.32) 



(2.33) 
(2.34) 



The C's are specific heats that would be obtained for the respective solutions. Since C a > C r , the free energies will 
intersect at a temperature T c _ < T\ (see Fig. 1), giving a first order transition at T c _. Ignoring the small differences 
between Tq, T c+ , and T± compared to Tq itself, we obtain 



AT C = T c+ - T c _ = 



fi — 1 



{Tc+ - TO, 



(2.35) 



where 



(i = (CJC r ) 1 ' 2 . (2.36) 
Differentiation of the free energies yields the entropies S r and S a , and the latent heat I = T(S r — S a )- 
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i = (jjL- l)C r AT c . 



(2.37) 



Let us ask if the scenario that has emerged could describe the actual data. To this end, we take AT C to be the 
observed T c splitting and C r and C a to be the specific heat jumps at the upper and lower transitions (relative to the 
normal state). Using data for samp le 1 from Ref. ||, AT C = 60 mK, C r — 98, and C a = 113 mJ/K mole UPt3, which 
implies fj, = 1.073. Equation (2.37) then gives I — 0.43 mJ/mole UPt3. If we conservatively take the experimental 
temperature resolution to be 10 /iK, then since the base specific heat is about 0.1 J/K mole, we conclude that a latent 
heat of 1 /iJ/mole would have been detected. The latent heat predicted by the 2D models is much larger. Further, 
the measured [g H c2 slopes of - 5 T/K imply values n = 6.6 x 10 3 KA 2 and V d = 73 mK. To obtain a AT C of 60 mK 
we would therefore need e ~ 70 mK, which is not much smaller than Vd- 



III. SIMPLE PERIODIC MODEL 



Tha analysis so far has been done assuming e -C Vd- It is also interesting to ask what happens when e>Kj. In 
this case we expect to obtain the local solution to a first approximation. Suppose the solution is of the strictly local 
type near T c+ . Consider two nearby domains with n||x, and ff||y, and an intervening domain with n at 120° to y 
and if _L n. The system will lower its gradient energy by tunneling of the r\ y component through the intervening 
domain. Since the ff component parallel to n is exponentially small in the ratio e/Vd, we expect the magnitude of the 
gradient energy gain will also be similarly small. However, the gradient energy favors relatively real r\ x and rj y , while 
the fa term favors a relative phase of tt/2 between them. General symmetry and continuity arguments would suggest 
that we still have two phase transitions with the lower one of first-order, but now the T c -splitting, and the degree of 
axiality at the lower transition will be exponentially small. It is then conceivable that an appreciably axiality only 
develops as a crossover at a lower temperature. This raises the possibility that the two true phase transitions are so 
close together that they appear to be one experimentally, and the crossover appears to be the phase transition at a 
lower temperature. 

To investigate the above possibility, we have numerically studied a toy model with a periodic antiferromagnetic 
domain structure. The assumption of periodicity eliminates the possibility of a first order transition altogether. This 
is because the model has an extra symmetry, and so is not in the same strict universality class as that studied in Sec. 
II. While the periodic model can not capture the two very close- by transitions (it only yields one), it can capture the 
low temperature crossover. Since we are only interested in the latter behaviour in this section, the periodic model 
is adequate. At the same time the periodicity simplifies the numerical problem enormously. We further take n(r) 
to vary only in one direction (x), and alternate between two orientations x and y with a periodicity of 2a. More 
concretely, we have 

^sbf= / dx2e(x)(\ Vx \ 2 ~\r, y \ 2 ), (3.1) 

where e(x) = ±e for < ±x < a, and e(x) — e(x + 2a). The restrictions of one-dimensonality and two orientations for 
n simplify the problem without changing its essential aspects. A closer analog of the modified 2D model || would be 
obtained by allowing k to vary periodically in addition to e, but since the fit to H C 2 data requires the two values of k 
to be rather close, the gradient energies are not appreciably different from domain to domain, and nothing is gained 
from this refinement. 

The eigenvalue problem for T c is identical to the Kronig-Penney model, and is analytically solvable. The solution 
is doubly degenerate corresponding to either r\ x or r\ y being zero everywhere. (Note that with Vd = n/a 2 as before, 
T c — Tq ~ e 2 /12Vd for small e, so that even with e ~ 3Vd, say, one is not in the strong symmetry breaking regime.) 
We solve the nonlinear GL equation below T c by discretizing it with a mesh size h > 2~ w a, and using a Newton- 
Raphson method. The free energy is computed for a closely spaced set of T values and numerically differentiated to 
obtain C v /T. The calculations are done with fa = 0.5Vd, and fa = 0.35Vd. This ratio of fa to fa is chosen to yield 
comparable C v /T jumps at T c ± in the uniform case. Our results for C v /T are shown in Fig. 2 for various e/Vd- There 
is clearly no second transition.^ For e/Vd — 500 there is a crossover that could at first sight mimic the data of Ref. 
||, but this is not so. As noted earlier, the measured || H c2 slopes imply Vd = 73 mK. Since the temperature axis 



8 This conclusion is supported by analytic work for e ^ Vd, in which T c is found by keeping only one wavevector q — n/a. A 
linear analysis about this solution below T c , and in the same approximation of one q, yields no further instabilities, except to 
the complex conjugate solution at T c itself. 
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in Fig. 2 has been scaled by e, and e = 36 K, the implied T c splitting would be 2.3e ~ 80 K. This is so far off the 
measured 50 mK splitting that the cross-over scenario within a 2D model can not be entertained for UPt3 . 

The surprising result from this exercise is that one might have expected a crossover and an apparent T c splitting of 
order e to appear for e = 20Vd, say. This is not so, and the crossover is only manifest at a much larger value of e/Vd- 



IV. DISCUSSION 



We have considered the effect of the short ranged antiferromagnetic domain structure on the superconducting 
transition in UPt3 within the 2D models, wherein the antiferromagnetism is the agency that splits the transition. It 
is known Jj"5| that when the AFM is weak and random, one does obtain two transitions but the lower one is of first 
order. We have shown that this prediction is in quantitative conflict with observations by calculating the latent heat 
at this transition, and showing that if the T c splitting and the specific heats are to agree with the observed data, 
then this latent heat is well above experimental limits on detectability. In the course of this analysis we have studied 
a variational axial state to better understand the short-ranged order. We have also studied the possibility that the 
lower transition may be a cross-over and not a true phase transition. In this case, the observed T c splitting can not 
be reconciled with the gradient energy scale. In either case, the AFM domain structure rules out all 2D models ^,|). 

Finally, we briefly disuss experiments which show that the two transitions in UPt3 merge into one at a pressure 
of p c ~ 4 kbar |fl9|| . The loose interpretation that this happens because the antiferromagnetism disappears at about 
the same pressure is thermodynamically unsustainable Q. All mean field models, whether based on one J^j or two 
J|,U primary order parameers, require either the reemergence of a splitting for p > p c , or a first order line at p ~ p c . 
Neither possibility has been ruled out so far in our view, as the broad heat capacity anomalies seen in Jl9| could 
conceal two transitions, and a first order transition has not been looked for. Careful experimental study of these 
possibilities is clearly highly desirable. At present, we believe that the two-order-parameter models, especially the 
AB model 0] , are the only viable ones left for describing UPt 3 . 
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FIG. 1. Sketch of the free energies (relative to that of the normal state) of the real and axial solutions, showing a change in 
slope, i.e., a first order transition at T c _. The transition at T c+ is of second order. 

FIG. 2. Computed specific heat for the simple periodic model, e and Vd are the symmetry breaking field and the gradient 
energies. 
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